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Abstract. Quantitative constraints on the wind clumping of massive stars can be obtained from 
the study of the hard X-ray variability of SFXTs. In these systems, a large fraction of the hard X- 
ray emission is emitted in the form of flares with typical duration of 3 ksec, frequency of 7 days and 
luminosity of 10 36 ergs/s. Such flares are most probably emitted by the interaction of a compact 
object orbiting at ~ 10 i?* with wind clumps (10 22-23 g). The density ratio between the clumps and 
the inter-clump medium is 10 2-4 . The parameters of the clumps and of the inter-clump medium 
are in good agreement with macro-clumping scenario and line-driven instability simulations. 



1. Introduction 

Stellar winds have profound implications for the evolution of massive stars, on the chemi- 
cal evolution of the Universe and as a source of energy and momentum in the interstellar 
medium. Indirect measures of the structure of massive-star winds are possible in X-ray bina- 
ries through the analysis of the interaction between the compact companion and the stellar 
wind. In this report we summarize the constraints obtained on wind clumping in HMXBs 
using t he hard X-ray variab ility observed by the IBIS/ISGRI instrument on board INTE- 
GRAL (jWinkler et al. 20031 ). Classical wind-fed, Roche-lobe underflow, super-giant HMXB 
(sgHMXB) are made of a compact object orbiting within a few (1.5 to 2.5) stellar radii from 
a super-giant companion. Recently INTEGRAL almost tripled the number of sgHMXB sys- 
tems known in the Galaxy and revealed a much more complex picture with two additional 
families of sources: (1) the highly-absorbed systems which have orbital and spin periods sim- 
i lar to those of clas sical sgHMXB but much higher absorbing column densities on average 



( IWalter et al. 20061 ) and (2) the fast transient systems which are c haracterized by fast ou t- 



bursts and by a very low quiescent luminosity (ISguera et al. 20061 : iNeeueruela et al. 20071 ). 



Several sources have now been proposed as candidate super-giant fast X-ray transient 
based on their hard X-ray variability characteristics, and, for a subset of them, optical 
counterpart spectral type. Contrasting statements have however been made on specific 
sources for what concerns their persistent or transient nature. In the frame of the current 
study (jl) we have considered all SFXT candidates together with several persistent and 
absorbed super-giant HMXB for comparison. 

We analyzed the available INTEGRAL data for 12 candidate SFXT that have large 
variability factors and compared them with the classical and absorbed sgHMXB systems that 
have a typical variability factor ^ 20. The sources have been separated in two categories. The 
SFXT include systems featuring hard X-ray variability by a factor ^ 100. "Intermediate" 
systems are candidate SFXT with smaller variability factors that could be compared with 
those of classical systems. 
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2. Hard X-ray Flares and Clumpy Winds 

The average count rate observed during SFXT flares lies between 3 and 60 ct/s which 
translates to hard X-ray luminosities of (0.2 — 4) x 10 36 erg/s. Such luminosities are not 
exceptional for sgHMXB but very significantly larg er than the typical X- ray luminosity of 
single massive stars of 10 30-33 erg/s at soft X-rays ( Cassinelli et al. 198ll ). As the sources 



are flaring at most once per day, their average hard X-ray luminosity is very low. It is 
therefore very unlikely that those systems have an average orbital radius lower than 10 13 cm, 
i.e. ~ 10 R*. One expects orbital periods larger than 15 days and underflow Roche lobe 
systems (note that no orbital period has yet been derived in any of these systems). 

The interaction of a compact object with a dense clump formed in the wind of a massive 
companion leads to increased accretion rate and hard X-ray emission. The free-fall time 
from the accretion radius R a = 2 x 10 10 cm towards the compact object is of the order of 
(2 — 3) x 10 2 sec. The infall is mostly radial (down to the Compton radius) and proceeds at 
the Bondi-Hoyle accretion rate. With a duration of tfi = 2 — 10 ksec, the observed short hard 
X-ray flares are significantly longer than the free-fall time. The flare duration is therefore 
very probably linked with the thickness of the clumps which, for a clump radial velocity 
Vd = 10 8 cm/s, is h d = V d x t fl ~ (2 - 10) x 10 11 cm. 

The average hard X-ray luminosity resulting from an interaction between the compact 
object and the clump can be evaluated as Lx = e M acc c 2 /tfi (where e ~ 0.1) and the mass 
of a clump can then be estimated as M d = (R d /R a ) 2 M acc = (R d /R a ) 2 Lx tfi/(e c 2 ) where 
R c i is the radius of the clump perpendicular to the radial distance. In the case of a spherical 

clump, M cl = ( IP ^) (^) 3 7.5 x 10 21 g. 

If N is the rate of clumps emitted by the star, the observed hard X-ray flare rate is 
given by T _1 = N(R 2 J AR 2 orb ). The rate of mass-loss in the form of wind clumps can then 

be estimated as M d = (iM^-^) ( T |^) 2 3 x lO" 6 M /y. For a (3 = 1 velocity 
law and spherical clumps, the number of clumps located between 1.05.R* and R or b can be 
evaluated as iV = (jj^j (jtjr) ( lo^cm ) ^.8 x 10 3 . Assuming spherical clumps, the clump 
density at the orbital radius is p d = ( 10 36 L ^ rg / s ) 7 x 10~ 14 g cm -3 and the corresponding 
homogeneous wind density is p h = M d /{A-n R 2 orb V d ) = ( 10 T d 10 s^ rg/s 3 * 1.5 X 10~ 15 g cm" 3 . 
The clump volume filling factor at the orbital radius is fv — ^ — {^j^j^) 0-02 and the 
corresponding porosity length is h = j^- = f^iTd) 15 x 10 12 cm. If the density of a clump 
decreases with radius as r~ 2/3 and its mass remains constant, the averaged homogeneous 
wind density within R obs isp^= NM d /(lnR 3 orb ) = (^ i^^ ^) 7 x 10~ 15 g cm" 3 and 
the average clump volume filling factor and porosity length could be estimated as 0.1 and 
3 x 10 12 cm, respectively. 

The variety of tfi, T and Ffi that are observed probably reflects a range of clump 
parameters and orbital radii. Several of the average clump parameters estimated above, in 
particular the clump density, filling factor and porosity length do not depend on the orbital 
radius, which is unknown, and only slowly depend on the observed quantities. 



3. Discussion 



The average clump parameters derived above match the macro-clumping scenario of (p) 
to reconcile clumping and mass-loss rates. The number of clumps derived above is also 
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comparable to evaluations by (jU; H). The volume filling factor, porosity length and the 
clump mass-loss rate are also similar to those derived by (j2|) from the study of ultraviolet 
and optical line profiles in two super-giant stars. 

The variation of the observed X-ray flux between flares and quiescence provides in 
principle a direct measure of the density contrast between the wind clumps and the inter- 
clump medium. Density contrasts of > 10 2 ~ 4 and 15-50 have been observed in SFXT 
and "Intermediate" sources, respectively. The density contrast is larger in SFXT than in 
"Intermediate" and, of course, classical systems. Density contrasts are probably stronger 
when clumping is very effective. 

Numerical simulations of the line driven instability (0) predict density contrasts as large 
as 10 3-5 in the wind up to large radii. At a distance of 10 i?*, the simulated density can 
vary between 10~ 18 and 10~ 13 g cm -3 and the separation between the density peaks are of 
the order of i?*. These characteristics are comparable to the values we have derived. 

Classical sgHMXB are characterized by small orbital radii R or b = (1.5 — 2.5) i?*, and by 
flux variability of a factor ^ 10. Such variabilities were modeled in terms of wind inhomo- 
geneities largely triggered by the hydrodynamic a nd photo-ionisation e ffects of the accreting 
object on the companion and inner stellar wind ( Blondin et al. 199ll ; Q]). At small orbital 



radii, the companion is close to fill its Roche lobe, which triggers tidal streams. In addi- 
tion the X-ray source ionizes the wind acceleration zone, prevents wind acceleration and 
generates slower velocities, denser winds, larger accretion radius and finally larger X-ray 
luminosities. Whether or not the stellar wind is intrinsically clumpy at low radius, the effect 
of the compact object on the wind is expected to be important. 

The main difference between SFXT and classical sgHMXB could therefore be their 
orbital radius. At very low orbital radius (< 1.5 R*) tidal accretion will take place through 
an accretion disk and the system will soon evolve to a common envelope stage. At low orbital 
radius (~ 2 i?*) the wind will be perturbed in any case and efficient wind accretion will lead 
to copious and persistent X-ray emission (10 36 ~ 37 erg/s). At larger orbital radius (~ 10 R*) 
and if the wind is clumpy, the SFXT behavior is expected as described above. If the wind 
clumps do not form for any reason, the average accretion rate will remain too low and the 
sources will remain mostly undetected by the current hard X-ray survey instruments. 
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